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Heavy-light decays such as B — > niv, B -> K*^ and B — >■ K^*^££ can be used to constrain the 
parameters of the Standard Model and in indirect searches for new physics. While the precision of 
experimental results has improved over the last years this has still to be matched by equally precise 
theoretical predictions. The calculation of heavy-light form factors is currently carried out in lattice 
QCD. Due to its small Compton wavelength we discretize the heavy quark in an effective non- 
relativistic theory. By formulating the theory in a moving frame of reference discretization errors 
in the final state are reduced at large recoil. Over the last years the formalism has been improved 
and tested extensively. Systematic uncertainties are reduced by renormalizing the m(oving)NRQCD 
action and heavy-light decay operators. The theory difi'ers from QCD only for large loop momenta 
at the order of the lattice cutoff and the calculation can be carried out in perturbation theory as 
an expansion in the strong coupling constant. In this paper we calculate the one loop corrections 
to the heavy-light vector and tensor operator. Due to the complexity of the action the generation 
of lattice Feynman rules is automated and loop integrals are solved by the adaptive Monte Carlo 
integrator VEGAS. We discuss the infrared and ultraviolet divergences in the loop integrals both 
in the continuum and on the lattice. The light quarks are discretized in the ASQTad and highly 
improved staggered quark (HISQ) action; the formalism is easily extended to other quark actions. 

PACS numbers: 12.15.Hh, 12.38.Bx, 12.38.Gc, 12.39.Hg, 13.20.He 



I. INTRODUCTION 

Decays of mesons containing heavy quarks provide an 
excellent laboratory for studying the heavy flavor sector 
of the Standard Model. For inclusive decays techniques 
such as quark-hadron duality and the operator product 
expansion have been used to predict decay amplitudes 
and spectra to high precision (see, for example Refs. 
[l|, m ) . Exclusive decays have a well defined hadronic fi- 
nal state and are easier to measure but obtaining precise 
theoretical predictions for these processes is more chal- 
lenging since the quarks in the final state are bound inside 
a hadron. Nevertheless improving on these predictions is 
crucial to check results from inclusive measurements and 
overconstrain the parameters of the Standard Model to 
uncover the effects of putative new physics. 

The decay B iriv can, for instance, be used to con- 
strain the magnitude of Vub, one of the least known 
Cabibbo-Kobayashi-Maskawa (CKM) mixing matrix el- 
ements. In addition, rare decays like B — >■ K*j or 
B — >• K^*Hi are loop-suppressed in the Standard Model 
and expected to be sensitive to the presence of new 
physics 
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Precise calculations of hadronic matrix elements are 
needed to match experimental precision (which is still 
to improve further with new results from LHCb and the 
planned SuperB factory) . Lattice QCD provides a model- 
independent framework for these calculations. It is, how- 
ever, difficult to discretize the heavy valence quarks di- 
rectly on lattices that are currently available since their 
Compton wavelength is comparable to the lattice spac- 
ing, a. To overcome this problem we use the effective lat- 
tice field theory NRQCD ^-Q, which describes QCD in 
the non-relativistic limit where the heavy-quark velocity 
is much less than unity. NRQCD is useful for studying 
mesons with heavy-quark constituents such as J/V'j 
D and B. The heavy-light vector form factor has been 
calculated in lattice NRQCD and combined with exper- 
imental data to extract \Vub\ [131 • "^^^ hadronic form 
factors are functions of the squared momentum transfer 
= (Pb — Pf)^ , where pB and pp are the momenta of 
the decaying B meson and the hadronic final state. Cur- 
rent lattice calculations using NRQCD work well only for 
large q^, partially owing to large discretization errors in 
the hadronic final state. The radiative decay B — > K*j 
has = 0, however, and most experimental data for 
B — > nil' comes from the small region. The large ex- 
trapolation to small q^ is a sizable source of systematic 
error in analyzing such decays and it is thus desirable 
to extend the range of accessible q^ in a lattice calcula- 
tion. In our approach, this is achieved by formulating the 
theory in a moving reference frame [ll| . The frame veloc- 
ity is chosen to reduce the three-momentum of the final 
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hadron state at small q , and hence also suppress the 
associated lattice artifacts. This approach is known as 
moving NRQCD (mNRQCD) ^^m,- The momentum of 
the heavy quark p = mw+prcs is split into a contribution 
mu associated with the frame velocity u = (7,77;) and 
a residual momentum prcs, which is of the order of the 
hadronic scale Aqcd- The first contribution is treated 
exactly so that corrections can be expanded in powers 
of Pres/"i- Discretization errors, which scale with some 
power of aprosi can be removed systematically by improv- 
ing the action. 

Although we focus here on mNRQCD, wc note that 
similar calculations arc possible using other descriptions 
of moving heavy quarks [igI - IiH . A comparison of pre- 
dictions from these different approaches would provide 
a valuable test of our understanding of the systematic 
errors in all these methods. 

Hadronic matrix elements of heavy-light currents have 
been calculated using lattice mNRQCD by the HPQCD 
Collaboration 0, [20|. The NRQCD effective theory 
is obtained from QCD by integrating out the effects of 
physics on energy scales of order m or larger; on the 
lattice this upper energy scale is determined by the lat- 
tice spacing, a. Operators in QCD are written as a for- 
mal expansion of operators defined in terms of the effec- 
tive non-relativistic fields of NRQCD. The coefficients 
of the NRQCD operators in the expansion are deter- 
mined as power series in the strong coupling constant 
as{q*/a) which is defined in an appropriate renormal- 
ization scheme and where q* is a dimension parameter 
of order unity. The NRQCD operators are ordered in 
powers of 1/m. For quarkonium matrix elements this 
gives a series in as and the quark relative velocity w^^^j 
whereas for heavy-light systems the series is an expansion 
in as and Aqcd /™- These coefficients are the radiative 
corrections which compensate for the missing contribu- 
tion of ultraviolet modes in QCD which are omitted by 
the imposition of the high-momentum ultraviolet cutoff 
in the lattice formulation. QCD and the effective lat- 
tice theory agree for infrared scales sufficiently below the 
heavy quark mass m and so the radiative corrections are 
governed by momenta larger than m > l/o where the 
strong coupling constant is small. It is therefore legiti- 
mate to calculate these corrections in perturbation the- 
ory. The coefficients are computed by equating the ma- 
trix elements of the operator in QCD and of its NRQCD 
expansion for any choice of external states. Because the 
coefficients are independent of the matrix element cho- 
sen, we carry out the radiative matching calculation for 
appropriately chosen external on-shell quark states; this 
is the usual technique employed in matching calculations 
for NRQCD (see, for example, [H [H). These radia- 
tive corrections are expected to be comparable in size to 
higher order corrections in the 1/m expansion. 

Owing to the complexity of the lattice actions used, the 
generation of Feynman rules has been automated |23l - [25j 
and the resulting integrals are calculated using the adap- 
tive Monte Carlo integrator VEGAS possibly after 



smoothing the integrand by adding an infrared subtrac- 
tion function. In jl5| radiative corrections to the heavy 
quark action have been calculated for a range of frame 
velocities. In this work we extend this calculation to the 
matching coefficients for the leading order operators of 
the vector and tensor currents. 

The outline of this paper is as follows. We discuss 
the different heavy-light continuum operators that con- 
tribute to i? — > Tr£iy and rare B decays in Sec. |TT1 and we 
introduce the various quark actions used in this project 
in Sec. IIIIl The central part of this work, the matching 
calculation between continuum and lattice operators, is 
presented in Sec. IIVI where we calculate one-loop matrix 
elements both in the continuum and in the effective lat- 
tice theory and combine them to obtain the matching 
coefficients. In Sec. IVlwe present numerical results. We 
summarize and discuss our findings in Sec. IVII 

Some preliminary results from this work appeared in 
Ref. filMlli. 



II. CONTINUUM OPERATORS 

(T) 

In general we shall denote current operators by J,\ ' 
with the appropriate number of Lorentz indices and 
where F labels the symmetry and transformation prop- 
erties and the integer subscript, n, distinguishes different 
operators with the properties. 

An example is the semileptonic decay B — > irliy which 
occurs at tree level in the quark picture. It is mediated 
by the hadronic vector current 



J. 



(1) 



where q and ^ are the light- and heavy-quark spinors. 
We follow the convention in Appendix A of Ref. |15| for 
the Dirac gamma matrices 7^. Only left-handed parti- 
cles (denoted by the subscript L) participate in the weak 
interaction. For massless leptons the total width of this 
decay is proportional to the square of the form factor 
f+{q^) defined by 



Ml -Ml 
q' ' 



(2) 



and the CKM matrix element Vub- 

The situation is more complicated for rare B decays 
which can only occur at loop level in the Standard Model. 
After integrating out physics at the electroweak scale, the 
transition is described by a set of effective operators Qj 
with their associated Wilson coefficients Cj{fi) [1^. For 
the radiative decay B — >■ K*'y the Hamiltonian is 



n 



(3) 
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Ci 


0.016 


C5 


0.017 


C2 


0.711 


Cs 


0.009 




-0.078 


C7 


-0.300 


C4 


0.093 


Cs 


-0.144 



TABLE I. Numerical values of the Standard Model Wilson 
coefficients in the effective Lagrangian ([3]) for the radiative 
decay B — > K*^. The coefficients Cj{fib) are evaluated 
in the leading logarithmic (LL) approximation at the scale 
jib = 5.0 MeV using the so-called "magic numbers" in Ref. [2^ 
with mz = 91.1876 MeV, mw = fJ.w = 80.425 MeV and 
as{mz) = 0.118. 



where the operator basis O, used in this work is given, 
for example, in Refs. (sol. Isij: 



(^1 = 

Q2 = 

Q3 = 

Qi = 



Q8 



{cLpTbLa) {sLal^iCLp ) , 
{cLaj'^bLa){SL/3^p.CLI3), 

q—u.d,c,s.b 

q—u,d,c,s,b 

q—u.d,c,s.b 

q—u.d,c^s,b 

e 



(4) 
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mbSLaCF'^''bRaF^u, 

2 mbSLaCr^''T^i3bRi3G1j^^, 



with 



(5) 



(G^j,) is the electromagnetic (chromodynamic) field 
strength tensor; a and (3 are color indices. 

This factorization separates the physics at large energy 
scales, which is contained in the model-dependent Wilson 
coefficients Cj (n) , from the universal hadronic matrix el- 
ements of the effective operators Qj. 




FIG. 1. Local tensor operator diagram contributing to the 
decay B — > K*^ 



of Cj are now known at next-to-next-to-leading (NNLL) 



order in the Standard Model [ij 



The Wilson coefficients are model dependent. For ex- 
ample, in [32| it is shown how C'/{fiw) and CsifJ-w) 
change in a two Higgs doublet model. The numerical 
size of these changes depends on the parameters of the 
specific model. In [s^ it is reported that the inclusive 
decay rate B — >■ X^j, which in the LL approximation is 
proportional to ICj {fj,b)\'^ , could in principle be enhanced 
by about a factor of three compared to the Standard 
Model. 

Although this enhancement is now ruled out by recent 
calculations of the Standard Model branching ratio in 
the inclusive B — > Xsj [H and exclusive B K*^ 
(33 - l36| decays, which are compatible with experimental 
results [131, the situation is less clear for the time de- 
pendent CP asymmetry in the exclusive decay: although 
it is expected to be small in the Standard Model [s^ it 
has not yet been measured to sufficiently high precision 
jsol l40j . In the Standard Model the opposite chirality 
operator, which is obtained by replacing 6_r i-> 6l and 
sl ^ SR\n Qj is suppressed by a factor of nislrrib as the 
spin fiip requires the insertion of a mass term. This is 
not necessarily the case in the new physics models stud- 
ied in [4l| where it is shown that the opposite chirality 
operator and mixing induced CP asymmetries can be en- 
hanced even if the branching ratio agrees with Standard 
Model predictions. 

Similar conclusions can be drawn for the decay 
B — > K'-*Hi [11,111], the forward-backward asymmetry is 
dependent on C7, Cg and Ciq. The current experimen- 
tal measurements [44l - l46| will be improved by the LHCb 
experiment [47j and help constrain these coefficients. 



Wilson coefficients 



By matching the effective Hamiltonian in Eqn. ([3|) to 
the Standard Model it is easy to see that at tree level 
C2 = 1 and all other coefficients are zero. After sum- 
ming the leading logarithms of the form [as \og{^w I fJ-b)] 
the dominant contribution to C7 comes from the mixing 
of Q2 to Qy via one loop diagrams [11]. The Standard 
Model Wilson coefficients relevant for the radiative de- 
cay are given in Tab. U in the leading logarithmic (LL) 
approximation. As can be seen from these numbers, the 
dominant operators are C2 and C7. Numerical values 



B. Local and non-local operators 

The operators in Eqns. pi4p can be split into two 
groups: the four quark operators 1-6 couple two hadronic 
currents at one point in space time. All other operators 
couple a heavy-light quark current to a gauge boson or a 
Icptonic current. These two sets of operators contribute 
differently to hadronic heavy-light matrix elements [48j . 




FIG. 2. Non-local four quark operator diagram contributing 
to the decay B — !> K*^ 



1. Local contributions 

The local contribution to the radiative decay B K*j 
is described by the tensor current (see Fig. [Ij 



167r2 



(6) 



The local contributions to the leptonic decay B — > K^* 
also include operators 



Qio 



I 



(7) 



which couple the heavy-light vector current J^^^ to a 
vector (or axial vector) leptonic current. 

In the following we will consider the vector current in 
Eqn. ([T} and the heavy-light hadronic tensor current 



J, 







(8) 



The heavy quark mass is included in this current as only 
left-handed particles participate in the weak interaction: 
flipping the chirality on one of the external legs requires 
the insertion of a mass term TO^2,9i? + (h.c.). Whilst in the 
Standard Model the operator mg(g^(T''''\E'L) with oppo- 
site chirality is suppressed by mq/m relative to Eqn. (|8|). 
this is not necessarily the case in new physics models |41| . 

As we set the light quark mass to zero in the match- 
ing calculation we will drop all chiral projectors in the 
following. 



2. Non-local contributions 

Non-local contributions come from diagrams like the 
one in Fig. [2] The gauge bosons couple to an internal 
quark loop which is created by contracting the two charm 
fields of a four quark operator. Given the size of C2 it 
is important to estimate the effect of these diagrams on 
the hadronic matrix element. 



FIG. 3. Resonant contribution to 6 
operators. 



S7 from four quark 



Long distance effects. The dominant contribution to 
the long distance amplitude induced by the h sec op- 
erators is usually assumed to come from the diagram 
where the photon couples directly to the charm quark 
loop. This is confirmed by the perturbative calculation 
combined with a quark model in [ioj where it is found 
that the main contribution generated by Q2 is given by 
diagrams where the photon couples to the cc loop and 
the gluon connects this loop to either the & or s quark. 
Only contributes at tree level but in [i^ it is argued 
that the 0{as) contribution of the four quark operator 
is of the same order as the one loop correction to the 
electromagnetic tensor operator Q^. 

The resonant contribution from the charm loop can 
be described as the decay B — )- Vtpn, where the ipn is a 
bound cc vector state, such as J/ip, which subsequently 
decays into a photon, see Fig. |31 

In this approximation the long-distance amplitude can 
be written as 



A 



{e;,)^{o\crcmA{B ^ ^ 



M2 



where e„ is the polarization of the vector meson ipn with 
mass M„ and width r„. For real photons with q^ = 
the sum is dominated by the lowest lying resonances. The 
mass of the J/ip is 3.097 GeV so that long distance effects 
from charm loops are expected to be suppressed by the 
inverse of this mass. 

This argument is supported by the explicit calculation 
in [sO] where the charm quark is integrated out to ob- 
tain an effective bsgj operator suppressed by The 
matrix element of this local operator is calculated us- 
ing QCD sum rules and found to be small, contributing 
around 5% of the dominant amplitude from Qy. 

Chromomagnetie tensor operator. The contribution 
from the chromomagnetie tensor operator is estimated 
There the decay amplitude is calculated for 
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both the electromagnetic and the tensor operator in the 
framework of a quark model. The contribution of Qg 
is found to be suppressed relative to Qy by a factor 
AqcD/n^B X Cs/Cj w 5%. 

To summarize, matrix elements of four quark oper- 
ators and of the chromomagnetie tensor operator are 
suppressed for small q^ and it is likely that they con- 
tribute little to the radiative decay. With currently avail- 
able techniques these non-local contributions can not be 
treated in lattice QCD and, if required, must be cal- 
culated using different approaches, such as QCD sum 
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rules [35|, |49| - |54| to complement the lattice calculation 
of the local operators. In the following we will concen- 
trate on the vector current Qj^j^^'^l and the tensor cur- 
rent ■m{qj^a^^'^'i> pj) which are the hadronic parts of the 
electroweak operators in ([6]) and Nonperturbative 
lattice matrix elements elements of these currents have 
been calculated [l^, H^l by using mNRQCD as an effec- 
tive theory for the heavy quark. 



III. LATTICE QUARK ACTIONS 

On currently available lattices the Compton wave- 
length of the b quark is smaller than the lattice spac- 
ing. We discretize the heavy quark in an effective, non- 
rclativistic theory where the high frequency fluctuations 
have been inte grat ed out. The construction of this action 
is described in [l5| and in the following we summarize the 
main results. 

The light quarks are discretized using highly improved 
relativistic actions [HI, [5^ . In the one loop calculations 
presented in this paper vacuum polarization effects of the 
light quark do not contribute. 



A. Moving NRQCD 

The (tree level) NRQCD action is obtained by decou- 
pling the quark and antiquark degrees of freedom in the 
fermionic action by a Foldy-Wouthuysen-Tani transfor- 
mation. The theory can be Lorentz-transformed to a 
moving frame and higher order time derivatives are re- 
moved by a subsequent field transformation. On the lat- 
tice, where Lorentz invariance is broken by discretiza- 
tion, this will give rise to a new theory which is known as 
m(oving) NRQCD. As the theory contains only first or- 
der time derivatives, propagators can be computed very 
efhciently by a single sweep through the lattice. 

The mNRQCD action is given by 

S^^^l {x, t) [iI;{x, t) - K{T)tlj{x, r - 1)] (10) 



operators defined in Ref. [l5|. In momentum space the 
non-relativistic dispersion relation 



with kernel 



The lowest order kinetic term Hq is 



Hc) = -iv ■ A 



2jm 



(12) 



E = v- p^.^ 



Pros - {V-Prcsf 



« i/m^ -I- {mjv + Pres)^ - 7"^ 
= \/m^ + p2 — 7m 



(13) 



is obtained. 5H contains higher order corrections in 1/m 
and operators which remove discretization errors. The 
action used in this work is described in Ref. [HI; it is 
correct to 0{l/m^, ^rd)' where froi is the relative velocity 
of the two quarks in a heavy-heavy system. 

The integer stability parameter n is introduced to re- 
move numerical instabilities for smaller quark masses m. 



B. Relativistic quark actions 

We separately consider two different staggered lattice 
actions describing the light quarks. 

The ASQTad action [55| suppresses "taste-breaking" 
interactions of lattice doublers by ©(aga^). This is done 
by introducing form factors for one gluon emission in the 
action. The Highly Improved Staggered Quark (HISQ) 
action reduces discretization errors further by an addi- 
tional level of smearing followed by link unitarization [sgI ] . 

Physically, the light quark mass is much smaller than 
the hadronic scale and in the matching calculation it will 
be set to zero. This simplifies the calculation and leads 
to additional relations between different matching coeffi- 
cients due to chiral symmetry. 

The MILC and UKQCD Collaborations have pro- 
duced a set of lattice configuration ensembles including 
ASQTad vacuum polarization effects and arc currently 
extending this to configurations with dynamical HISQ 
fermions [s^, HI] . These configurations are used in the 
nonperturbative calculation of heavy-light form factors 



IV. MATCHING CALCULATION 

To match heavy-light operators, expressed in terms of 
effective mNRQCD heavy quark fields, to those in the 
continuum theory, we must calculate radiative correc- 
tions both in the continuum and on the lattice. We match 
the theories at leading order in the 1/m expansion and 
one loop order in a^. Matrix elements of 0{l/m) oper- 
ators are expected to be of the same size as the leading 
radiative corrections and are matched at tree level. 



Continuum calculation 



where v is the frame velocity and A^, A^^) and aI^-* are In the continuum we calculate the one loop matrix 
first and second order gauge covariant finite difference elements of the currents Jq and expand the result in 
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inverse powers of the heavy quark mass. We use the no- 
tation T = V,T to denote the expressions in Eqns. (jllSp . 
Owing to Lorentz invariance these expansions can be ex- 
pressed as hnear combinations of a smah number of tree 
level matrix elements: 

At leading order in the 1 /m expansion the operators that 
contribute to the vector current are: 



J, 



(V)m 



and to the tensor current: 



(15) 



(16) 



where = (7,71?) is the frame velocity. The 

one loop contributions to the mixing matrix 

^jr.con) _ J otgbz'^J'"'^"^^ -t- . . . are for the vector 
operator [21| 



r^(V,con) 1 11 3 \2 / 2 



(y,con) 



37r' 



(17) 



and for the tensor operator 



(Txon) 




1 

3^ 



27 



■ log }? j-rr? + \ log 771^ / [T 



(18) 



We introduced a gluon mass A to regulate infrared di- 
vergences. It is important to use the same infrared reg- 
ulator in both continuum QCD and the effective lattice 
theory, any dependence on the gluon mass will cancel in 
the matching coefficients. 

As the tensor current is not conserved, its anomalous 
dimension does not vanish: 



(con) Sofs 



(19) 



The coefficients of the infrared logarithms in Eqns. P7)) 
and (jl8p agree due to heavy quark symmetry. 



B. Construction of lattice operators 

On the lattice we must construct operators Jq'"''^*-' 
which have the same on-shell matrix elements as the as- 
sociated continuum operators: 



At tree level, the operators in the effective theory are 
obtained from Eqns. (jllSp by applying the field transfor- 
mation 

*(x) = S'(A)f(5)e-™'^°"-^Ai5, — (21) 

5(A) is a spinorial Lorentz boost, T the FWT trans- 
formation decoupling the quark- and antiquark fields 
in the rest frame and Ad^ an additional field trans- 
formation to remove higher order time derivatives, 
^'^(x) = (V'«(x),0)"^ is the (positive energy) field in the 
effective theory. 

Using the explicit expressions in Ref. [l5l | one finds at 
©(l/m) 



/7 



170!) ■ T> ■ D iv ■ D 



2m 



2m 



27m 



5(A)*, 



(22) 

from which the tree level currents in the effective theory 
can be read off by inserting the field transformation in 
Eqns. (ITT81) . 

In the following we calculate the one loop matching 
coefficients of the leading order operators in the 1 /m ex- 
pansion. At this order both T and A^^ are equal to the 
identity. The Lorentz boost is 



5(A) = 



1 



^2(1 + 7) 



((1 + 7) - ^vi} ■ 77") . 



(23) 



We choose v, the direction of the frame velocity, to be 
along one of the lattice axes, i.e. v — (1,0,0). We 
can then classify the lattice directions (and associated 
Lorentz indices) as timelike (denoted 0, as usual), paral- 
lel to V (denoted ||) or perpendicular to v (denoted _L). 



1 . Choice of operator basis 



In the continuum the operator basis Jq , Ji is used. 
On the lattice Lorentz invariance is broken and it is con- 
venient to work in another basis which is spanned by op- 
erators with different Dirac structure. Firstly, the opera- 
tors Jq and are split into the sum of two operators 



(r) 



Forn^O: 



4? = P*'^Vi(")g(a:)rM'„, 

Jf} = -P^^^f2{v)q{x)rv ■ 77°*, . (24) 

In the vector case F = 7'^ and p^^^ = 1, and in the tensor 
case F = cr^'^ and p'^^-' ~ m. The velocity dependence 
has been absorbed in the functions 



hiv) 



1+7 . / ^ / 7 

J2[V) = V 



^7 



2(1+7) 



(25) 



The corresponding leading order (in the heavy quark 
(20) expansion) operators Jj are obtained by replacing 
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r 1-^ u^^ in the vector case, and F M- 2i{'y^u'^ — j^u^^) 
for the tensor. 

The tree level matrix elements of all these operators 
form a basis for expanding higher order matrix elements. 

(T) 

In the following we do not, however, need to include ^ , 
because the tree level matrix elements are given by linear 
combinations of those of Jq y due to heavy quark sym- 
metry ~ 'i'v In the vector case: 



jiV)0 

•^1,1 

7(^^)11 
•^1,1 



(1 



-7)4!^^" 



j{V}0 
''1,2 

7(^0 1! 

•^1,2 



7'^0,2 ' 



iv)\\ 



and for the tensor operator: 



t(T)0|1 
•^14 

7(T)0|| 
"^1,2 

t(T)0± 
■^1,1 

.(T)OJ- 
■^1,2 

7(r)||-L 
•^1,1 

7(T)||± 
"^1,2 



(l"7)4i 



(26) 



(7J™I + (1 ^ ,;.0,2 y 
((1-7)J™'-7J™'), 



-2(1 
-2(1 



7)^r" 



2 



(27) 



Clearly, this decomposition is not Lorentz invariant, but 
can be carried out for fixed frame velocity. 

On the lattice the two operators in Eqn. mix under 
renormalization. 



toZ^""'^™) and Zj""'^™): 



^^(V,co„)0 
^^(y,con)0 
^^(V,co„)|| 
^^(V,con)|| 
^^(y,con)± 



,5Z 



(V.con) 







(V.con) 



{V,con) 







szl 



(V^,con) 



-iszi 



(y,con) 



SZ 



(y,con)_ 



SZ, 

= 



(y,con) 



(31) 



For the tensor operator there is no dependence on the 
Lorentz indices: 



SZ 



(T,con) _ ^^(T,con) 



(32) 



2. Matching coefficients 



Combining Eqns. 



the lattice operator which 



has the same one loop matrix elements as the continuum 
operator is 



J, 



(r,lat) 



(l + a,cf ) Jo,+ 



with 



(r,lat) 



(33) 



(34) 



((?|j(^,)|6)trcc, (28) 

(with j,k = 1,2). Instead of using this basis (which we 
will call the (1, 2) basis), it is more convenient to work in 
the (+, — ) basis: 



7(r) _ 7(r) , j{r) 



(29) 



as only Jq _l contributes to processes at tree level. We 
may then write 

{q\4':l\b),on = il + aJZ^^-'°''^){q\4'^l\b)trcc+ 
a^SZ 



(30) 



(r.lat) 



in this ba- 



On the lattice renormalization factors SZ 
sis can then be defined in an analogous way to Eqn. 
[i.e. where j,k = ±]. For the vector operator, we must 
distinguish whether the Lorentz index of the current is 
timelikc, parallel or perpendicular to the frame veloc- 
ity. Relations like Eqn. can then be used to relate 



C. Mixing matrix 

In the (1, 2) basis of operators, the lattice mixing ma- 
trix can be split into a diagonal part and a contribution 
ifk from one particle irreducible (IPI) diagrams. 



SZ 



(r,lat) 

jk 



SZ. 



(r) 

mult 



SZf 



f(r) 



(35) 



For the vector current the multiplicative renor- 
malization contains the wavefunction renormalization 
only: SZ'^.^^^ = ^{SZq + SZ^). For the tensor cur- 
rent, however, there is an additional contribution 
from the renormalization of the heavy quark mass: 
5Z^2xt = \i.5Z, + SZ^,)-SZ,^. 

The relation between renormalized and bare param- 
eters is vr = ZyV, 7fl = (1 - t>|)"^/^, mR =^ Z„im 
and q = ^JWqqR, = ^JZ^,^r. mji is the pole mass, 
which can be defined perturbatively both in the contin- 
uum and on the lattice. The renormalization constants 
can be expanded perturbatively using the generic formula 
Zx = 1 + asSZx + 

The renormalization of the velocity functions /1.2 in 
Eqn. ([251) is f],R = Zf- fj with 



SZf 



^-—^SZ,, SZf, = ^^SZ,. (36) 
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One then finds 

.^(r.iat) _ . (r) 1 . (D 

^^(r>t) ^ 2jz^^^(r)^ (37) 

Even though we use only the leading order heavy-light 
operators we still include 1/m corrections in the action. 
Next we isolate infrared divergences in the renormaliza- 
tion constants and find (in Feynman gauge) 

SZq = ^loga^X^ +Fg, 

OTT 

2 

SZj, = -— loga^A^ + Fj,{v,am), 
on 

SZm = F,niv,am), 
SZj; = Fy{v, am), 

d+ = -^^og a'X' + F^ll {v, am), 
^^+l=F^l[iv,am). (38) 
I 



Here we make the lattice spacing a explicit and A is 
the gluon mass used in the continuum calculation. The 

infrared divergence of C++ is independent of the Dirac 
structure F due to heavy quark symmetry and can 
be inferred from the subtraction integral discussed in 
Sec. lIVHTl The functions Fx are infrared finite and can, 
if required, be expanded in powers of the inverse heavy 
quark mass on the lattice. 



It is interesting to note that the logarithms in Eqn. ((38)) 
represent both infrared (A — )- 0) and ultraviolet (a 0) 
divergences. A similar combination of short and long dis- 
tance divergences occurs in HQET if the theory is reg- 
ulated in dimensional regularization: The heavy quark 
propagator does not contain any scales and the integral 
vanishes in this case (see, for example, Ref. [1^). In the 
HQET case, this simplifies the calculation of the match- 
ing coefficients as only the QCD integrals need to be cal- 
culated. 



D. Results for the vector operator 



After cancelling infrared divergences the final expression for the matching coefficients is 

-\iF, + F„) + \Fy + ^ log a'm' F^^f , 
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^+ — 


^12^ 




1 


- 


47r 




11 




^12^ 




27 




37r 


,(_V)I1 ^ 


27 








2 



I ,^ „ . 1 „ 1 



2_2 Tp(V)|| 



-2^-9^ - '4' J -r ^Fv + — log a'm' - F^^^ , 
+ Fi>) + + ^ log - ' 

2 " «+- ' 

-If.-f^:i\ 

~F^2"- (39) 

(V) 

In the limit w — > the operator Jq 2 does not contribute as it is proportional to v; in NRQCD there is only one 
operator with matching coefficient c^'"-' : 

cf^ = c'+^ + cL^^ -+ c(^) for v^O. (40) 

We find 

c(^)" = -±-liF, + F„) + l- loga'm' ~ F^', 



An 2'" 2n 
c(^)^ = - i (^^, + F^) + ^ loga^m^ - F^P\ (41) 



The matching coefficient of the zero component of the vector (or axial vector) current at w = has been calculated 
in [gO] and the corresponding calculation for the spatial components can be found in [l^]. We tested our code by 
reproducing these results. 
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E. Results for the tensor operator 

The corresponding results for the tensor operator are: 



and in the NRQCD hmit u = 



F. The anomalous dimension 

The ultraviolet behavior of the lattice theory is de- 
scribed by the logarithmic terms in Eqns. (|39l42p . In 
particular, the log a^m^ term is a UV divergence which 
is independent of the Dirac structure of the rcnormalizcd 
operator due to heavy quark symmetry. As the short 
distance behavior of the effective theory is different from 
that of continuum QCD, its coefficient is not the same as 
that of the \ogm^ /jj? term in Eqn. The anomalous 

dimension of the lattice operator can be obtained by not- 
ing that the renormalized operator is related to the bare 
operator by multiplication by Zj 



(lat). 



^(r.rcn) _ /7(lat)\ ^ JT) 



(44) 



The counterterm has to be chosen such that it absorbs 



the logarithmic UV divergence in SZ 



(r,lat) 

-+ ' 



,(lat) 



1 - ^ [log"Vilt + (finite terms)] 



(45) 



where fii^t is an arbitrary scale which cancels in physical 
results. We thus find 



(lat) 

Tr 



1 dZ- 



(lat) 



as 

TT 



(46) 



This agrees with the result for HQET regularized in di- 
mensional regularization 61 1. 



G. Quark renormalization parameters 

The wavefunction renormalization of massless ASQTad 
quarks has been calculated to one loop in We re- 

peated this calculation with a larger number of points in 
the VEGAS integration to obtain 

-0.92411(42)+ -!- log a^A^ (47) 

in Feynman gauge (the error is the statistical error of the 
VEGAS integral). For massless HISQ quarks we obtain 



g ^ASQTad 



^7HiSQ ^ -0.3905(16) + — loga^A^ 



3tt 



(48) 
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Q 




FIG. 4. One particle irreducible diagram at one loop 



The one loop renormalization of the heavy quark wave- 
function, frame velocity and mass is given for a range of 
frame velocities in Ref. [l5|. As demonstrated there the 
magnitude of all renormalization parameters is reduced 
significantly by including mean fields corrections. Only 
the wavefunction renormalization has a logarithmic IR 
divergence, given by — -^loga^A^. 



H. One particle irreducible matrix elements 

The one particle irreducible (IPX) matrix elements can 
be found by evaluating the one loop diagram in Fig. |4l 
The 0{as) one particle irreducible correction to the 

operators Jq J = fj'qTSj'i'i, (with Si = 1, S2 = —7 • '£'70) 
is given by 

(<?l4j^|fc)iat,iPi = asfjuip')^Pu 



fe=l,2 

= a. E ^fk'fku{p')^SkU (49) 

fe=l,2 

where the heavy-quark four spinor is 

u">-{^:). w.,. .-.{(;),(?)}^ (SO) 

(V) 

is the lattice integrand, after factoring out external 
quark spinors, the strong coupling constant and the ve- 
locity function fj. To extract Qk j replace the spinors 
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by Euclidean [63 on-shell projection operators: 

u^nt= J2 u^"^u^-~>^ = \{l + 7o) (51) 



and take the trace of 



fc=l,2 



The Dirac matrix II'-'"-' is a suitable projection operator 
which depends on F; as both sides of (|52|) are a (linear) 
function of the four momentum p', this relation defines 
for anj,fc. 



1. Infrared subtraction function 

In some configurations in momentum space the inte- 
gration contour in the fco plane is pinched by poles. This 
leads to large peaks in the integrand and can generate 
infrared divergences in the final result. 



We construct an appropriate infrared subtraction func- 
tion ^(^'"ti) s]-|2ooth the integrand and thus speed up the 
convergence of the VEGAS estimate of the integral. The 
IPI matrix elements can be written as 



AT) 



(4" 

j.(sub) 



lat) 



^'(sub) 
Jjk 



(2^ 



-r- 



(sub) 



(53) 



Construction of the subtraction function is guided by the 
continuum integral, which has the same infrared behavior 
as the corresponding lattice expression. In the continuum 
the IPI correction to the operator Jg ■ at the matching 
(52) point p' = 0, p = (m, 0) is given by 



iljl 



Df\k){~gT'^v,) 



(54) 



where D^j^"^ (fc) is the heavy quark propagator at 
p = (m, 0). This integral can be rendered UV finite with- 
out changing the infrared structure by replacing 



<)(fc) 



kn — iv ■ k 



1-^ 



27TO 



(fc 



(55) 



where the metric is (7^,^ — diag(+l, +1, +1, +1) in Eu- 
clidean space, and the velocity four-vector is given by 
u = (17,71;). As in Eqn. (|49l) we write 



(r,sub) 



167r 



d^k 



E 



(2,)4"b')(-^7p) 

fk&:''^u{p')^s,u. 



"fc2" 



VS,\J 



{k+pf 



F + A2 



(56) 



The subtraction integral is independent of the Dirac 
structure and is easy to solve analytically: 

^(sub) _ f d'^k (sub) 

^^'^ ~ J (27r)4 ■'i'' 

= - (1 + log AVm2) + 0(A/m). (57) 

This concludes our discussion of the structure of the 
matching calculation. In the next section we present 
numerical values for different quark masses and frame 
velocities. 



V. NUMERICAL RESULTS 

In Figs. [5] and ini we show results for the matching co- 
efficients for the vector and tensor current (see Tables [ll] 



to |V] for numerical values) . 

In both cases we use a heavy mass of m = 2.8 and a 
stability parameter n = 2. These are the values currently 
used in nonperturbative calculations of heavy-light form 
factors on coarse MILC lattices d^H^. The gluon action 
is Symanzik improved [6^. We present results both for 
the ASQTad and HISQ light quark action. 

For the vector current we calculate the matching coef- 
ficients for three different directions of the Lorentz index 
/i: temporal (^ = 0), parallel to the frame velocity (de- 
noted II) and perpendicular (_L). The frame velocity is 
chosen to be along the lattice axis /i = 1. For w = we 
consider a corresponding set of directions: fi ~ 0, fi = 1 
and jjL = 2. 

For the tensor current there arc four different cases for 
indices {^i,v): (Q,||), (0,_L), (||,_L) and (-L,-L). For u 
we choose {^JL,v) = (0,1), (0,2), (1,2) and (2,3). The 
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0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



frame velocity v frame velocity v 




0.2 0.4 0.6 0.8 1 
frame velocity v 



FIG. 5. Matching coefficients, vector current (top) and ten- 
sor current (bottom). The heavy quark mass is m = 2.8 and 
the ASQTad action is used to discretize the hght quark. Due 
to chiral symmetry discussed in the main text the following 
matching coefficients are identical for the tensor current: c'}'' 
(open red squares) and cL^ (open red downward triangles), 
(open red circles) and ci^ (open red upward triangles), 
c!|J' (filled blue squares) and c^''^ (filled blue downward tri- 
angles), c?^ (filled blue circles) and c![.^ (filled blue upward 
triangles) . 

renormalization scale of the tensor current is fi = m. 



A. Discussion 

As the light quark is massless its propagator and vertex 
functions anticommutc with 7'^. With cr°^ = 17'^ cr^'^, this 
implies that matching coefficients for (0, ||) and (-L,_L) 
are identical (given that we boost in the direction fi = 1). 
The same holds for (0, _L) and (|j,-L), as cr°^ = —ij^a^^. 
For V = the matching coefficients for all combina- 
tions of {n, v) agree as there is no preferred direction. 
Heavy quark symmetry Iff'^'^v = ^v) and relations such 
as 7°7-' = —ia^^ can be used to relate IPI matrix ele- 
ments of the vector and tensor current. 

We emphasize that the magnitude of all matching co- 
efficients is reduced by including mean field corrections 
in the renormalization parameters, the dependence on 
the frame velocity is weak and all matching coefficients 




0.2 0.4 0.6 0.8 1 
frame velocity v 



FIG. 6. Matching coefficients, vector current (top) and tensor 
current (bottom). The heavy quark mass is m = 2.8 and the 
HISQ action is used to discretize the light quark. See the 
comment on Fig. [5] for relations between the tensor current 
matching coefficients. 



(except for c\_ , which depends on the renormalization 
scale) are of order 0.3 or smaller for moderate v. With 
as{2/a) sa 0.3 the relative size of the radiative corrections 
does not exceed 10%. From heavy quark power counting 
we expect the matrix elements of the 0(1/to) operators 
(which are matched at tree level) to be suppressed by 
the same factor Aqcd/"t- « 0.1 relative to the leading 
operators. 

The matching coefficients appear to diverge for u — 1. 
It should be noted that whereas mNRQCD reduces to 
NRQCD in the limit w — ?> 0, the coUinear theory at w = 1 
is qualitatively different, in particular both mass and ve- 
locity renormalization are not defined in this limit. 

In addition, for large v the momentum distribution of 
the heavy quark in the initial state meson is boosted by a 
factor of 7 so that for highly relativistic frame velocities 
the power counting in 1/m will break down. As shown 
in [l5| the statistical errors of nonperturbative matrix 
elements grow with decreasing and in practise it is 
unlikely that they will be calculated for v > 0.5. 
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1. Vector current 

The matching coefficient for the zero component of the 
vector current at m = 2.8, c^^)" ^ 0.04293(52) in Ta- 
ble ini is in perfect agreement with the corresponding 
value po = 0.043(2) in Table III of For u = wc 

find that the matching coefhcients c^^^^ = 0.26970(40) 
and c(^)2 ^ 0.26929(39) agree within errors as would be 
expected from rotational invariance and are consistent 
with = 0.270(1) in Table II of [H. 

The splitting between the matching coefficients for dif- 
ferent Lorentz indices is reduced by using the HISQ light 
quark action. This reduction is not, however, as pro- 
nounced as for the tensor current. In the continuum the 
matching coefficient c^Y^ is zero. Using the HISQ action 

reduces the magnitude of c^^" by nearly a factor of two 

and an even stronger reduction is observed for cf^^". 



2. Tensor current 

For the tensor current we find that the splitting be- 
tween the matching coefficients for (/i, i') = (0, ||) and 
(0, _L) is reduced by using the HISQ action for the light 
quark. A similar reduction is seen in the splitting be- 
tween the (II, -L) and (-L,-L) matching coefficients. 

The matching coefficients c^^^"^ are always very small 

and their magnitude is about 0.1. The size of c^^'"' 
depends on the continuum renormalization scale /i; for 
/i = m, we find that these coefficients are also very small 
when using ASQTad light quarks. They are, however, 
larger when the HISQ action is used to discretize the 
light quark. 



3. Quark mass dependence 

We repeated the calculations for a heavy quark mass 
m = 1.9, which corresponds roughly to the bare quark 
mass on the fine MILC lattices. The results for a range 
of frame velocities are shown in Tables IVII to IIXI The 
absolute size of the matching coefficients is larger than 
for TO = 2.8 but still typically lies in the range 0.1 - 0.5 
for the frame velocities considered. 



VI. CONCLUSION 

M(oving) NRQCD is a useful tool for extending the 
range of accessible in a lattice calculation of heavy- 
light form factors. The formalism has been improved 
and tested extensively over the last years. Radiative cor- 
rections to the effective actions have been calculated in 
a previous publication [isj . Further reduction of system- 
atic uncertainties is justified by an increase in precision of 



experimental results. In this paper we show how system- 
atic errors due to radiative corrections can be reduced 
by renormalizing the heavy-light vector and tensor cur- 
rents. After cancelling infrared divergences the one loop 
corrections to matching coefficients are of the order one 
and smaller. The results will be used in the current cal- 
culation of nonperturbative form factors [20j . 

As the lattice imposes a cutoff ^ 1/a, operators which 
are formally suppressed by 1/m can "mix down" to the 
leading order operators |6J] . These power law terms can 
be suppressed by constructing perturbatively subtracted 
0(1/to) operators jj^^^^'" = jj^^ _ agf^^ Jq^^ which do 
not mix down to the leading operators at one loop order. 

(V) 

f^.g is calculated from the one particle irreducible correc- 
tions to the 0(1/to) operators. We find that a.t v — 0.4 
the non-perturbative matrix element of the subtracted 
1/to operators is a factor of around 0.05 smaller than the 
leading order matrix element, which is consistent with 
heavy quark power counting where AgcD/'Ti ^ 0.1. Be- 
fore subtraction the ratio of the matrix elements can be 
as large as 0.3. 

In this work we used the ASQTad and HISQ actions 
to discretize the light quark but it should be noted that 
our approach is easily extended to other discretizations 
such as the Domain Wall fermion action [g^ . 

With currently available techniques only matrix ele- 
ments of local operators in rare exclusive decays can be 
computed in lattice QCD. Often these operators describe 
the dominant effects. For the radiative decay B — >■ K*j 
the four quark operators Q3, • . . , are suppressed by 
their small Wilson coefficients. At the physical point 
where = 0, the contribution of Q2 is suppressed as 
the cc vector resonance which connects this operator and 
the external photon is far off shell. Model calculations 
show that the matrix elements of the chromomagnetic 
operator Qs arc small. This implies that the dominant 
contribution comes from the electromagnetic tensor op- 
erator Qy and the matrix elements of this operator can 
be evaluated in lattice QCD. 

An additional complication is that the effective heavy 
quark theory is only valid at maximum recoil, i.e. at 

= m^. Results have to be extrapolated to = 
using a phenomenological ansatz. A simple, physically 
motivated parametrization is given in (66l . [67| . Recently 
a model independent parametrization, which uses the an- 
alyticity of the form factors, has been suggested in josj . 

Even if the dominant contribution to a given process 
is not given by a local operator, lattice calculations are 
still useful when combined with other approaches such as 
QCD sum rules. 
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Appendix A: Pole shift 

In this appendix we discuss the choice of contour for 
the lattice integrals. As for the self-energy calculations in 
Ref. [l^ , care has to be taken when choosing the integra- 
tion contour in the fco plane (where k is the momentum 
of the gluon in the loop). The heavy quark pole must lie 
inside the integration contour (in the z = e**^" plane) but 
for certain values of the loop momentum it lies outside 
the unit circle. The integration contour then has to be 
deformed to ensure that the result can be Wick-rotated 
back to Minkowski space. 

In the following we discuss the corresponding contour 
shift for the lattice three point functions. We begin by 
mapping the positions of the poles of the light quark 
propagator and then discuss the contour choice for the 
full integrands. 



1. Poles of the light quark action 

We use the ASQTad and HISQ actions to describe the 
light, relativistic quarks and, as discussed in the main 
text, we treat these quarks as being massless. The prop- 
agators for these actions are identical, with denominator 



A = ^sin^(fc^) (l + ^sm^{k^ 
= w(l + iw)2 + s(fc). 



(Al) 



where .s(fc) = ^^^j^ sin^(fcj) (l -f i sin^(fcj)) and 
u) = sin^(fco). 

The poles of the propagator correspond to A 0. 
For given fixed, spatial three-momentum k (with 



kj G [— TT, tt]), finding the poles reduces to solving a cubic 
equation with real coefficients, 

- I2y + 36.s(fc) - 16 = 0, (A2) 

where y = a; + 4. This equation either has three real 
solutions, or one real solution and one conjugate pair of 
complex solutions, depending on the sign of the discrim- 
inant D [6^: defined by 



£» = g2 - 64 with q = 18s(fc) - 8 



(A3) 



It has three real solutions if I? < (or cquivalently 
s{k) < 8/9): 



2/1 = -2P cos /3, y2,3 - 2P cos ± - j (A4) 

where f3 = ^ arccos((//P'^) and P = 2sgn(7 = ±2. 

Alternatively, it has one real solution and one conju- 
gate pair of complex solutions if D > (or cquivalently 
s{k) > 8/9): 

2/1 = -2Pcosh/3, 2/2,3 = P {cosh 13 ±iV3smh (3) (A5) 

where /? = ^ arccosh(g/P'^) and P is defined as above. 
In either case, for each a; = 2/ — 4, there are four solu- 



tions for 



, which can be labelled as: 



z±± = ± V 1 - 2a; ± 2v/u;2 _ 



(A6) 



Of these twelve poles, two are "physical" and survive in 
the continuum limit. They can be identified as those that 
lie on the unit circle |z| = 1 for fc ^ 0. The other ten 
"spurious" poles are lattice artifacts with masses propor- 
tional to and therefore decouple in the continuum 
limit. These poles (sometimes called ghost poles) are, of 
course, not the lattice doublers. 

We calculated z for a large number of random momenta 
fc with kj G [— TT, tt] ; the resulting distribution of the poles 
in the complex plane is shown in Fig. [T] 

It is useful to compare the positions of these poles of 
the improved actions with those of the naive propagator. 
The equation describing the position of the poles is linear 
in cu, leading to four solutions for z: 



(naive 

4± 




±2yk (i + fe ) 



(A7) 



with k = X]')=i ^i'^^l^)- show the comparison as 



a function of |fe| in Fig. [5] For each spatial momentum, 
all the poles z in the improved propagator lie outside 

(naive) 

4- 



the region defined by ' 



< \z\ < Iz'^Z""']- This 
allows us to use the positions of the naive poles as a safe 
lower bound for the positions of the improved poles when 
choosing integration contours. 
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Re(z) 



-2 



|zl = 1 



-6 -4-2 2 
Re {z} 



FIG. 7. Poles of the massless ASQ/HISQ fermion propagator 
in the complex plane. 
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E: 1 



0.1 




1 



1*1 



FIG. 8. Absolute value of poles in the massless naive and 
ASQ/HISQ (improved) fermion propagator as a function of 

|fc| = sin^(fcj). The plot is generated by finding the 

poles for a large random sample of momenta kj G [— vr, tt]. 



To take care of causality one can employ the ie pre- 
scription, changing the denominator to 



A = ^sin2(fc^) ( 1 + isin2(fc^,) ) ^ ie 



6 



(A8) 



For the discussion of the effects of this, it is sufficient to 
concentrate on the three poles outside the unit circle with 
positive real part (the other poles are related to these by 
transformations z —> l/z,—z,—l/z). One of these is the 
physical pole, which for s(fe) = lies on the unit circle 
with a small negative imaginary part. In addition, there 
are two additional, spurious poles with larger real part: 
one with a negative (and small) imaginary part and one 
with a positive (and small) imaginary part. 

The movement of these poles as the momentum in- 
creases is shown in Fig.[9l As s(fc) gets larger, the physi- 
cal pole moves outwards. The spurious poles, meanwhile, 
move in opposite directions: one moves outwards just 
below the real axis and away from the physical pole; the 



lzl=l 




Im(z) 



s(k)=8/9 



FIG. 9. Movement of three of the twelve poles in the mass- 
less ASQ/HISQ propagator in the z complex plane (the other 
poles are related by the transformations z 1/z, —z, —1/z). 
The arrows show the motion as s(fe) increases. The physi- 
cal pole is shown with a solid (red) line. The spurious poles 
are shown using a dashed (blue) line. The unit circle is also 
shown (in green). The contour pinch occurs for s{k) = 8/9. 



other, meanwhile, moves inwards just above the real axis, 
towards the physical pole. 

When s{k) — 8/9, the physical pole touches (within 
distance 2e) one of the spurious poles and both, now 
being complex conjugates of each other, start to move 
away from the real axis in opposite directions. 

Having established the positions of the poles, in the 
next section we will use these to choose the contours in 
the lattice integration appropriately. 



2. Pole shift in one particle irreducible integrals 

In this section we discuss the poles of one particle ir- 
reducible three point integrals in Sec. IIVHI 

Let the position of the heavy quark pole in the z plane 
be denoted by z^ and the poles of the naive gluon prop- 
agator by z± such that |z_| < 1 < \z+\. The two poles of 

the naive light quark action are z^f ^ , whereas the six poles 

of the improved light quark action are located at z^'^ 
(and the corresponding positions with opposite sign), or- 
dered such that 



|zi^)3|<|zi^)|<|zi^)|<l< 



4\\ 



< z 



W I 



(A9) 

Note that, as discussed in the previous section, only one 
of the poles is physical. 

From the calculation of heavy quark renormalization 
parameters, it is known for naive (Wilson) glue that 
\zh\ < p^ . The poles of the Symanzik-improvcd 
gluon action lie outside the band defined by z_ < \z\ < 
2;+, so the same holds for improved gluons. We can there- 
fore concentrate on the relative positions of the poles of 
the heavy and (improved) light quark propagator. 

At high frame velocities and for certain choices of spa- 
tial momentum, it turns out that the heavy quark pole 
can cross poles of the light propagator outside the unit 
circle. Note, however, that as discussed in the main 
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/ ^ ■■" X. 
/ "X. 
/ X 
// \ 


/ \ 


( > 


\ / 

\v /y 



-7C -7t/2 7l/2 7C 

X 

FIG. 10. Absolute values of poles as a function of x with 
fc = (a;, 0,0). Poles of naive gluons are shown with a dashed 
(green) line; of improved fermions with a long dashed (blue) 
line; of the simple mNRQCD action (Ho only) with a dotted 
(magenta) line; and of the full mNRQCD action with a solid 
(red) line. The frame velocity is w = 0.95, the heavy quark 
mass m = 2.8 and the stability parameter n = 2. 



text, it is unlikely that very large frame velocities will 
be used in the evaluation of non-pcrturbative matrix ele- 
ments. Examples are shown in Fig. I10[ where we choose 
fc = (.T, 0, 0) with — TT < X < TT and discuss both a simple 
action with only and the full mNRQCD action. 

The crossings are seen for certain negative values of x, 
where < \zh\ < k+.2l- The problem gets worse if 

the full mNRQCD action is used. To be able to Wick- 
rotate back to Minkowski space in these cases, the con- 
tour needs to be deformed such that it encloses the heavy 
quark pole but not the light quark poles outside the unit 
circle. Suitable contours are shown on the left in Fig. [11] 
In a similar way to that used in Ref. [22|, the contours 
can be deformed to avoid the poles as much as possi- 
ble, arising at the triple contours shown on the right in 
Fig.ini 

Computationally, the procedure is as follows. We 
choose the contour(s) separately for each value of the 
spatial momentum (generated, for instance, by the VE- 
GAS integration code): 

1. \z}i\ < z_, As Z}i is not the smallest negative 
pole and the contour does not need to be shifted 
from \z\ = 1. 

2. < \zh\ < z+,|z^^\|. The contour 
is shifted outwards to halfway between \zh\ and 
min{z+,zy.^}. 

3. \zfj <\zn\< fori = 1, 2 (see Fig.HO]). 
A pole crossing has occurred and it is necessary to 
integrate along three contours: (a) anticlockwise 
without shift, \z\ = 1; (b) clockwise, shifting the 
contour midway between and \zh\] and (c) 



counterclockwise with the contour between \zji \ and 
min{|z^^_^^._^J,z+}. 

To speed up the vegas calculation, in cases (1) and (2) 
we first check using the poles of the naive light quark ac- 
tion, only calculating the poles of the improved version 
if that test is inconclusive. When shifting a contour mid- 
way between poles Za and Zf,, we shift the contour to 

As the pole crossing only occurs for large momenta this 
is a lattice artifact which would disappear in the contin- 
uum limit. It must, however, be included in a lattice- 
continuum matching calculation. 
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FIG. 11. Integration contour for momentum space configurations with pole crossing. Tfie light quark poles are denoted by x, 
the naive gluon poles by o and the heavy quark pole by □. 
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V c(^')° c(^')i c(^')2 

0.00 0.04293(52) 0.26970(40) 0.26929(39) 







^+ 


^+ 


,10 


0.1945(21) 


-0.0553(22) 


0.1384(18) 


,20 


0.1965(12) 


-0.0500(13) 


0.1411(12) 


,30 


0.2005(10) 


-0.0419(11) 


0.1444(10) 


,40 


0.20708(84) 


-0.0293(10) 


0.15071(88) 


,50 


0.21660(77) 


-0.01038(88) 


0.15972(84) 


,60 


0.22620(73) 


0.01358(85) 


0.16775(85) 


,70 


0.24297(71) 


0.05349(80) 


0.17817(87) 


,75 


0.25596(72) 


0.08412(80) 


0.18243(91) 


,80 


0.27674(77) 


0.12920(85) 


0.1862(10) 


,85 


0.3179(13) 


0.2057(13) 


0.1918(12) 


,90 


0.4244(15) 


0.3604(16) 


0.1986(15) 


,95 


0.6492(19) 


0.6615(19) 


0.2186(23) 





,(V)0 






,10 


-0.1533(22) 


0.3227(22) 


0.1301(17) 


,20 


-0.1547(13) 


0.3220(13) 


0.1274(10) 


,30 


-0.1598(10) 


0.3183(10) 


0.12224(81) 


,40 


-0.16838(86) 


0.31106(85) 


0.11325(72) 


,50 


-0.18137(81) 


0.29887(80) 


0.09967(67) 


,60 


-0.19460(80) 


0.28291(80) 


0.08532(68) 


,70 


-0.21033(81) 


0.25390(80) 


0.06649(69) 


,75 


-0.21513(84) 


0.23150(83) 


0.05808(73) 


,80 


-0.21139(91) 


0.20128(91) 


0.05445(80) 


,85 


-0.1931(11) 


0.1530(12) 


0.0561(10) 


,90 


-0.1353(13) 


0.0713(14) 


0.0777(12) 


,95 


0.0615(23) 


-0.0832(23) 


0.1347(22) 



TABLE II. Vector current matching coefficients for heavy quark mass m = 2.8. The ASQTad action is used to discretize the 
light quark. The Lorentz indices can be timelike (0), parallel (||) or perpendicular (1.) to the frame velocity. The table shows 
the central value and error from the VEGAS integration. 
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y f,^T)0,l ^{T)0,2 g(T)l,2 g(T)2,3 

0.00 0.0762(12) 0.0761(12) 0.0763(12) 0.0761(12) 





„(T)0,|| 


{T)0,X 






,10 


-0.0507(24) 


-0.0159(21) 


-0.0157(22) 


-0.0463(24) 


,20 


-0.0559(17) 


-0.0206(16) 


-0.0214(16) 


-0.0526(17) 


,30 


-0.0675(16) 


-0.0325(15) 


-0.0321(15) 


-0.0652(16) 


,40 


-0.0861(15) 


-0.0496(15) 


-0.0497(15) 


-0.0846(15) 


,50 


-0.1148(16) 


-0.0748(15) 


-0.0750(15) 


-0.1136(16) 


,60 


-0.1605(17) 


-0.1135(16) 


-0.1136(16) 


-0.1595(17) 


,70 


-0.2302(19) 


-0.1666(18) 


-0.1668(18) 


-0.2295(19) 


,75 


-0.2791(21) 


-0.1984(21) 


-0.1992(21) 


-0.2782(21) 


,80 


-0.3379(26) 


-0.2280(25) 


-0.2292(25) 


-0.3368(26) 


,85 


-0.4006(34) 


-0.2397(35) 


-0.2397(35) 


-0.3990(34) 


,90 


-0.4463(53) 


-0.1688(53) 


-0.1696(53) 


-0.4441(53) 


,95 


-0.279(14) 


0.224(14) 


0.223(14) 


-0.279(14) 





^(T)O.II 






(T)±,± 
C_ ' 


,10 


0.1184(20) 


0.0849(18) 


0.0871(18) 


0.1189(20) 


,20 


0.1147(11) 


0.0839(11) 


0.0855(11) 


0.1149(11) 


,30 


0.10781(85) 


0.07967(86) 


0.08054(86) 


0.10774(85) 


,40 


0.09607(74) 


0.07144(76) 


0.07214(77) 


0.09583(74) 


,50 


0.07857(69) 


0.05896(72) 


0.05912(72) 


0.07848(69) 


,60 


0.05889(69) 


0.04414(72) 


0.04412(72) 


0.05860(69) 


,70 


0.03355(70) 


0.02162(73) 


0.02142(73) 


0.03332(70) 


,75 


0.02163(74) 


0.00812(77) 


0.00787(76) 


0.02131(74) 


,80 


0.01460(82) 


-0.00494(84) 


-0.00522(84) 


0.01425(82) 


,85 


0.0146(10) 


-0.0199(10) 


-0.0205(10) 


0.0142(10) 


,90 


0.0361(13) 


-0.0310(13) 


-0.0312(13) 


0.0357(13) 


,95 


0.0980(22) 


-0.0125(22) 


-0.0114(22) 


0.0976(22) 



TABLE III. Tensor current matching coefficients for heavy quark mass m = 2.8. The ASQTad action is used to discretize the 
light quark. 
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V 




c(v)i 




0.00 


-0.10173(91) 


-0.03811(87) 


-0.03825(87) 




V 




^+ 




0.10 


-0.1323(22) 


-0.1980(22) 


-0.1593(20) 


0.20 


-0.1318(15) 


-0.1942(15) 


-0.1575(14) 


0.30 


-0.1286(13) 


-0.1902(13) 


-0.1513(12) 


0.40 


-0.1226(12) 


-0.1835(12) 


-0.1417(12) 


0.50 


-0.1147(11) 


-0.1732(12) 


-0.1276(11) 


0.60 


-0.1079(11) 


-0.1636(11) 


-0.1116(11) 


0.70 


-0.0982(11) 


-0.1454(11) 


-0.0895(11) 


0.75 


-0.0906(11) 


-0.1307(11) 


-0.0760(12) 


0.80 


-0.0775(11) 


-0.1067(11) 


-0.0610(12) 


0.85 


-0.0465(12) 


-0.0578(12) 


-0.0412(14) 


0.90 


0.0268(13) 


0.0411(14) 


-0.0157(16) 


0.95 


0.2011(22) 


0.2624(21) 


0.0229(23) 




V 








0.10 


0.0322(21) 


0.1574(20) 


0.1218(16) 


0.20 


0.0303(13) 


0.1562(12) 


0.1183(10) 


0.30 


0.0269(10) 


0.15322(93) 


0.11199(78) 


0.40 


0.02062(84) 


0.14692(82) 


0.10099(70) 


0.50 


0.01122(78) 


0.13632(76) 


0.08488(66) 


0.60 


0.00253(77) 


0.12358(76) 


0.06641(67) 


0.70 


-0.00754(78) 


0.10023(77) 


0.04171(68) 


0.75 


-0.00961(81) 


0.08269(80) 


0.02950(72) 


0.80 


-0.00407(88) 


0.05947(88) 


0.02188(79) 


0.85 


0.0132(10) 


0.0209(10) 


0.01954(93) 


0.90 


0.0648(13) 


-0.0436(13) 


0.0357(12) 


0.95 


0.2373(22) 


-0.1771(23) 


0.0933(22) 



TABLE IV. Vector current matching coefficients for heavy quark mass m — 2.8. The HISQ action is used to discretize the light 
quark. 
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y f,^T)0,l ^{T)0,2 f.mi.2 g(T)2,3 

0.00 -0.2317(15) -0.2318(15) -0.2316(15) -0.2317(15) 





„(T)0,|| 


{T)0,X 






,10 


-0.3610(25) 


-0.3330(23) 


-0.3305(23) 


-0.3555(25) 


,20 


-0.3643(19) 


-0.3380(18) 


-0.3375(18) 


-0.3614(19) 


,30 


-0.3741(17) 


-0.3500(17) 


-0.3496(17) 


-0.3726(17) 


,40 


-0.3890(17) 


-0.3691(17) 


-0.3690(17) 


-0.3881(17) 


,50 


-0.4120(17) 


-0.3974(17) 


-0.3970(17) 


-0.4109(17) 


,60 


-0.4480(18) 


-0.4410(18) 


-0.4414(18) 


-0.4469(18) 


,70 


-0.5022(21) 


-0.5040(20) 


-0.5043(20) 


-0.5011(21) 


,75 


-0.5394(23) 


-0.5437(22) 


-0.5441(22) 


-0.5384(23) 


,80 


-0.5829(27) 


-0.5843(26) 


-0.5854(26) 


-0.5827(27) 


,85 


-0.6264(35) 


-0.6102(34) 


-0.6103(34) 


-0.6264(35) 


,90 


-0.6497(53) 


-0.5754(53) 


-0.5761(53) 


-0.6490(53) 


,95 


-0.458(14) 


-0.240(14) 


-0.240(14) 


-0.459(14) 





^mo.ii 






(T)±,± 
CI 


,10 


0.1230(19) 


0.0964(17) 


0.0955(18) 


0.1234(19) 


,20 


0.1187(11) 


0.0939(11) 


0.0938(11) 


0.1187(11) 


,30 


0.11138(82) 


0.09018(84) 


0.09064(84) 


0.11126(82) 


,40 


0.09914(72) 


0.08398(75) 


0.08432(75) 


0.09898(72) 


,50 


0.08099(67) 


0.07380(70) 


0.07417(70) 


0.08075(67) 


,60 


0.05968(67) 


0.06302(70) 


0.06342(70) 


0.05937(67) 


,70 


0.03107(68) 


0.04637(71) 


0.04668(71) 


0.03060(68) 


,75 


0.01635(73) 


0.03647(75) 


0.03669(75) 


0.01583(72) 


,80 


0.00591(80) 


0.02745(82) 


0.02778(82) 


0.00529(80) 


,85 


0.00098(94) 


0.0168(10) 


0.0169(10) 


0.00026(94) 


,90 


0.0146(12) 


0.0102(12) 


0.0108(12) 


0.0141(12) 


,95 


0.0705(22) 


0.0299(22) 


0.0302(22) 


0.0696(22) 



TABLE V. Tensor current matching coefficients for heavy quark mass m = 2.8. The HISQ action is used to discretize the hght 
quark. 
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V 


c{v)o 






0.00 


-0.06534(59) 


0.33741(41) 


0.33681(41) 




V 




^+ 


^+ 


0.01 


0.246(18) 


-0.185(18) 


0.162(14) 


0.10 


0.2164(21) 


-0.1954(22) 


0.1514(18) 


0.20 


0.2147(13) 


-0.1924(14) 


0.1501(12) 


0.30 


0.2170(10) 


-0.1849(12) 


0.1506(10) 




V 




c(_v)ll 




0.01 


-0.294(18) 


0.530(18) 


0.185(14) 


0.10 


-0.2831(21) 


0.5337(22) 


0.1848(17) 


0.20 


-0.2829(13) 


0.5353(13) 


0.1847(11) 


0.30 


-0.2889(11) 


0.5333(11) 


0.18050(89) 



TABLE VI. Vector current matching coefficients for heavy quark mass m — 1.9. The ASQTad action is used to discretize the 
hght quark. 



V 


c(T)0,l 






c(r)2,3 


0.00 


0.3410(16) 


0.3411(16) 


0.3414(16) 


0.3412(16) 




V 


„(T)0,ll 


_(T)0,-L 
^+ 


^+ 


^+ 


0.01 
0.10 
0.20 
0.30 


0.141(18) 
0.1524(26) 
0.1448(21) 
0.1319(19) 


0.210(14) 
0.2094(24) 
0.2017(20) 
0.1907(18) 


0.205(14) 
0.2085(24) 
0.2017(20) 
0.1903(18) 


0.167(18) 
0.1577(26) 
0.1476(21) 
0.1339(19) 




V 


(r)o,|| 
c_ 




g(T)||,± 




0.01 
0.10 
0.20 
0.30 


0.162(18) 
0.1784(20) 
0.1777(12) 
0.1712(10) 


0.119(14) 
0.1268(18) 
0.1267(12) 
0.12228(94) 


0.119(14) 
0.1285(18) 
0.1273(12) 
0.12256(94) 


0.190(18) 
0.1814(20) 
0.1785(12) 
0.17183(95) 



TABLE VIL Tensor current matching coefficients for heavy quark mass m = 1.9. The ASQTad action is used to discretize the 
hght quark. 
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V 


c{v)o 






0.00 


-0.15556(94) 


0.01365(88) 


0.01296(88) 




V 




^+ 


^+ 


0.01 


-0.095(17) 


-0.273(17) 


-0.164(13) 


0.10 


-0.1180(23) 


-0.2782(22) 


-0.1721(20) 


0.20 


-0.1198(16) 


-0.2779(16) 


-0.1720(15) 


0.30 


-0.1176(13) 


-0.2744(14) 


-0.1689(13) 




V 




c(_v)ll 


,(_v-)x 


0.01 


-0.055(17) 


0.283(17) 


0.189(13) 


0.10 


-0.0359(21) 


0.2908(21) 


0.1840(17) 


0.20 


-0.0366(13) 


0.2926(13) 


0.1831(11) 


0.30 


-0.0415(10) 


0.2914(10) 


0.17766(88) 



TABLE VIII. Vector current matching coefficients for heavy quark mass m = 1.9. The HISQ action is used to discretize the 
light quark. 



V 


c(T)0,l 




c(T')1.2 


c(r)2,3 


0.00 


0.0171(17) 


0.0175(17) 


0.0172(17) 


0.0173(17) 




V 


„(T)0,ll 


^+ 


^+ 


^+ 


0.01 


-0.199(17) 


-0.123(13) 


-0.112(13) 


-0.174(17) 


0.10 


-0.1800(27) 


-0.1163(25) 


-0.1150(25) 


-0.1763(27) 


0.20 


-0.1863(22) 


-0.1231(21) 


-0.1229(21) 


-0.1840(22) 


0.30 


-0.1966(20) 


-0.1347(20) 


-0.1345(20) 


-0.1950(20) 




V 


(r)o,|| 
c_ 




g(T)||,± 




0.01 


0.208(17) 


0.135(13) 


0.125(14) 


0.183(17) 


0.10 


0.1915(19) 


0.1285(17) 


0.1287(18) 


0.1889(19) 


0.20 


0.1889(12) 


0.1282(12) 


0.1287(12) 


0.1876(12) 


0.30 


0.18199(93) 


0.12518(94) 


0.12555(94) 


0.18114(93) 



TABLE IX. Tensor current matching coefficients for heavy quark mass am = 1.9. The HISQ action is used to discretize the 
light quark. 
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